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Abstract. The linear anisotropies in the temperature of the cosmic microwave 
background (CMB) radiation and its polarization provide a clean picture of fluctua- 
tions in the universe some 370 kyr after the big bang. Simple physics connects these 
fluctuations with those present in the ultra-high-energy universe, and this makes 
the CMB anisotropies a powerful tool for constraining the fundamental physics that 
was responsible for the generation of structure. Late-time effects also leave their 
mark, making the CMB temperature and polarization useful probes of dark energy 
and the astrophysics of reionization. In this review we discuss the simple physics 
that processes primordial perturbations into the linear temperature and polariza- 
tion anisotropies. We also describe the role of the CMB in constraining cosmological 
parameters, and review some of the highlights of the science extracted from recent 
observations and the implications of this for fundamental physics. 



1 Introduction 

The cosmic microwave background (CMB) radiation has played an essential 
role in shaping our current understanding of the large-scale properties of the 
universe. The discovery of this radiation in 1965 by Penzias and Wilson pp, 
and its subsequent interpretation as the relic radiation from a hot, dense phase 
of the universe [2] put the hot big bang model on a firm observational footing. 
The prediction of angular variations in the temperature of the radiation, due 
to the propagation of photons through an inhomogeneous universe, followed 
shortly after but it was not until 1992 that these were finally detected 
by the Differential Microwave Radiometers (DMR) experiment on the Cos- 
mic Background Explorer (COBE) satellite The fractional temperature 
anisotropies are at the level of 10 -5 , consistent with structure formation in 
cold dark matter (CDM) models [51 E], but much smaller than earlier predic- 
tions for baryon-dominated universes PJEJ. Another experiment on COBE, 
the Far InfraRed Absolute Spectrophotometer (FIRAS), spectacularly con- 
firmed the black-body spectrum of the CMB and determined the (isotropic) 
temperature to be 2.725 K 
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In the period since COBE, many experiments have mapped the CMB 
anisotropics on a range of angular scales from degrees to arcminutes (see ^U] 
for a recent review) , culminating in the first-year release of all-sky data from 
the Wilkinson Microwave Anisotropy Probe (WMAP) satellite in February 
2003 \H] ■ The observed modulation in the amplitude of the anisotropies with 
angular scale is fully consistent with predictions based on coherent, acous- 
tic oscillations 7 , derived from gravitational instability of initially adiabatic 
density perturbations in a universe with nearly-flat spatial sections. The am- 
plitude and scale of these acoustic features has allowed many of the key cos- 
mological parameters to be determined with unprecedented precision |12| , and 
a strong concordance with other cosmological probes has emerged. 

In this review we describe the essential physics of the temperature anisotropies 
of the CMB, and its recently-detected polarization ^2]> and discuss how these 
are used to constrain cosmological models. For reviews that are similar in 
spirit, but from the pre- WMAP era see e.g. |14U15| . We begin in Sect.|21with 
the fundamentals of CMB physics, presenting the kinetic theory of the CMB in 
an inhomogeneous universe, and the various physical mechanisms that process 
initial fluctuations in the distribution of matter and spacetime geometry into 
temperature anisotropies. Section|21discusses the effect of cosmological param- 
eters on the power spectrum of the temperature anisotropies, and the limits 
to parameter determination from the CMB alone. The physics of CMB polar- 
ization is reviewed in Sect. 01 and the additional information that polarization 
brings over temperature anisotropies alone is considered. Finally, in Scct.|3we 
describe some of the scientific highlights that have emerged from recent CMB 
observations, including the detection of CMB polarization, implications for in- 
flation, and the direct signature of dark-energy through correlations between 
the large-scale anisotropies and tracers of the mass distribution in the local 
universe. Throughout, we illustrate our discussion with computations based on 
ACDM cosmologies, with baryon density fl^h 2 = 0.023 and cold dark matter 
density fl c h 2 = 0.111. For flat models we take the dark-energy density pa- 
rameter to be Oa = 0.75 giving a Hubble parameter H = 73kms _1 Mpc -1 . 
We adopt units with c = 1 throughout, and use a spacetime metric signature 
+ • 

2 Fundamentals of CMB Physics 

In this section we aim to give a reasonably self-contained review of the essential 
elements of CMB physics. 

2.1 Thermal History and Recombination 

The high temperature of the early universe maintained a low equilibrium frac- 
tion of neutral atoms, and a correspondingly high number density of free elec- 
trons. Coulomb scattering between the ions and electrons kept them in local 
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kinetic equilibrium, and Thomson scattering of photons tended to maintain 
the isotropy of the CMB in the baryon rest frame. As the universe expanded 
and cooled, the dominant element hydrogen started to recombine when the 
temperature fell below ~ 4000 K - a factor of 40 lower than might be an- 
ticipated from the 13.6-eV inoization potential of hydrogen, due to the large 
ratio of the number of photons to baryons. The details of recombination are 
complicated since the processes that give rise to net recombination occur too 
slowly to maintain chemical equilibrium between the electrons, protons and 
atoms during the later stages of recombination |161 IT7| (see JS| I0r recent re- 
finements) . The most important quantity for CMB anisotropy formation is the 
visibility function - the probability that a photon last scattered as a function 
of time. The visibility function peaks around ~ 370 kyr after the big bang, 
and has a width ~ 115 kyr, a small fraction of the current age ~ 13.5 Gyr |12j . 
After recombination, photons travelled mostly unimpeded through the in- 
homogeneous universe, imprinting fluctuations in the radiation temperature, 
the gravitational potentials, and the bulk velocity of the radiation where they 
last scattered, as the temperature anisotropies that we observe today. A small 
fraction of CMB photons (current results from CMB polarization measure- 
ments ^HI indicate around 20 per cent; see also Sect. 15. if underwent further 
scattering once the universe reionized due to to the ionizing flux from the first 
non-linear structures. 

2.2 Statistics of CMB Anisotropies 

The spectrum of the CMB brightness along any direction n is very nearly 
thermal with a temperature T(n). The temperature depends only weakly on 
direction, with fluctuations AT(n) at the level of 10 -5 of the average tem- 
perature T = 2.725 K. It is convenient to expand the temperature fluctuation 
in spherical harmonics, 

AT(n)/T = ai m Yi m (n) , (1) 

irn 

with a\ m = (— l) m a;_ m since the temperature is a real field. The sum in 
Q runs over I > 1, but the dipole (I = 1) is usually removed explicitly 
when analysing data since it depends linearly on the velocity of the observer. 
Multipoles at / encode spatial information with characteristic angular scale 

The statistical properties of the fluctuations in a perturbed cosmology 
can be expected to respect the symmetries of the background model. In the 
case of Robertson- Walker models, the rotational symmetry of the background 
ensures that the multipoles a\ m are uncorrelated for different values of I and 
m : 

{ai m a*, m ,) = Ci5u>5 mm r , (2) 

which defines the power spectrum C;. The angle brackets in this equation 
denote the average over an ensemble of realisations of the fluctuations. The 



4 



Anthony Challinor 




• • VSA 

• WMAP 

■ □ CBI Mosaic 
ACBAR 



Fig. 1. Compilation of CMB anisotropy measurements (as of February 2004) 
from WMAP (black filled circles), the Very Small Array (VSA gD]; shaded cir- 
cles representing two interleaving binning schemes) , the Cosmic Background Imager 
(CBI 1211 122| ; open and filled squares for two different binning schemes) and the 
Arcminute Cosmology Bolometer Array Receiver (ACBAR |23|: triangles). (Figure 
reproduced, with permission, from |20|.) 



simplest models of inflation predict that the fluctuations should also be Gaus- 
sian at early times, and this is preserved by linear evolution of the small fluc- 
tuations. If Gaussian, the ai m s are also independent, and the power spectrum 
provides the complete statistical description of the temperature anisotropies. 
For this reason, measuring the anisotropy power spectrum has, so far, been 
the main goal of observational CMB research. Temperature anisotropies have 
now been detected up to I of a few thousand; a recent compilation of current 
data as of February 2004 is given in Fig. ^ 

The correlation between the temperature anisotropies along two directions 
evaluates to 

O] i 1 

(AT(ni)AT(n 2 )) = T 2 ^ -^—CiPi(cos9) , (3) 

i 

which depends only on the angular separation 9 as required by rotational 
invariance. Here, Pi(x) are the Legendre polynomials. The mean-square tem- 
perature anisotropy is 

(AT 2 ) = T 2 ^2 ~r~~Ci « T 2 f l Stp± Cl A\nl , (4) 
i 71 J k 

so that the quantity 1(1 + l)C;/27r, which is conventionally plotted, is approx- 
imately the power per decade in I of the temperature anisotropies. 



2.3 Kinetic Theory 



The CMB photons can be described by a one-particle distribution function 
f(x a ,p a ) that is a function of the spacetime position x a and four- momentum 
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p a of the photon. It is defined such that the number of photons contained 
in a proper three-volume element d 3 x and with three-momentum in d 3 p is 
/d 3 a;d 3 p. The phase-space volume element d 3 a;d 3 p is Lorentz-invariant and is 
conserved along the photon path through phase space (see, e.g. [23)- ^ follows 
that / is also frame-invariant, and is conserved in the absence of scattering. 
To calculate the anisotropies in the CMB temperature, we must evolve the 
photon distribution function in the perturbed universe. 

To avoid over-complicating our discussion, we shall only consider spatially- 
flat models here, and, for the moment, ignore the effects of polarization. For a 
more complete discussion, including these complications, see e.g. j2Sl|2S|- Cur- 
vature mostly affects the CMB through the geometrical projection of linear 
scales at last scattering to angular scales on the sky today, but has a negligible 
impact on pre-recombination physics and hence much of the discussion in this 
section. The subject of cosmological perturbation theory is rich in method- 
ology, but, for pedagogical reasons, we adopt here the most straightforward 
approach which is to work directly with the metric perturbations. This is also 
the most prevalent in the CMB literature. The 1+3-covariant approach is 
a well-developed alternative that is arguably more physically-transparent than 
metric-based techniques. It has also been applied extensively in the context 
of CMB physics H3 EH El EDI ED E2 • The majority of our discussion will 
be of scalar perturbations, where all perturbed three-tensors can be derived 
from the spatial derivatives of scalar functions, although we discuss tensor 
perturbations briefly in Sect. 12.51 

For scalar perturbations in spatially-flat models we can choose a gauge 
such that the spacetime metric is [35] 



where rj is conformal time (related to proper time t by dt = adrj), a is the scale 
factor in the background model and, now, x is comoving position. This gauge, 
known as the conformal Newtonian or longitudinal gauge, has the property 
that the congruence of worldlines with constant x have zero shear. The two 
scalar potentials <j> and ip constitute the scalar perturbation to the metric, 
with <p playing a similar role to the Newtonian gravitational potential. In the 
absence of anisotropic stress, <fi and ip are equal. We parameterise the photon 
four- momentum with its energy e/a and direction e (with e 2 = 1), as seen by 
an observer at constant x, so that 



Free photons move on the geodesies of the perturbed metric, p^V^p 1 ' = 0, so 
the energy and direction evolve as 



ds 2 = a 2 (ry)[(l + 2^)dr/ 2 - (1 - 2^)da; 2 ] , 



(5) 



= a" 2 e[l - V,(l + 0)e] . 



(6) 



de/di] 
de/dr/ 



edip/dij + e(<ft + ip) , 

v±(0 + v) , 



(7) 

(8) 
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where dots denote d/dr, and Vj_ is the three-gradient projected perpendicular 
to e. We see immediately that e is conserved in the absence of perturbations, 
so that the energy redshifts in proportion to the scale factor in the background 
model. The change in direction of the photon due to the projected gradient 
of the potentials in the perturbed universe gives rise to gravitational lensing 
(see e.g. for a review). 

The dominant scattering mechanism to affect CMB anisotropics is classical 
Thomson scattering off free electrons, since around recombination the average 
photon energy is small compared to the rest mass of the electron. Furthermore, 
the thermal distribution of electron velocities can be ignored due to the low 
temperature. The evolution of the photon distribution function in the presence 
of Thomson scattering is 

^ = _ (1 + VOw/ + t|-o(1 + </>Ra T / /(e, e')[(l + (e • e'f] de' 
dr, lo7r J 

df 

- an c a T e ■ v h e— , (9) 
oe 

where n B is the electron (proper) number density, or is the Thomson cross 
section, and the electron peculiar velocity is = dx/dr,. The derivative on 
the left of is along the photon path in phase space: 

t- = if + e ■ V ^ + ^ - e ' V ^) £ ?r ( 10 ) 

dr, or, oe 

to hrst order, where we have used J7J) and (JSJ and the fact that the anisotropies 
of / are first order. The first term on the right of JjJJ describes scattering out of 
the beam, and the second scattering into the beam. The final term arises from 
the out-scattering of the additional dipole moment in the distribution function 
seen by the electrons due to the Doppler effect. In the background model / is 
isotropic and the net scattering term vanishes, so that / is a function of the 
conserved e only: / = /(e). Thermal equilibrium ensures that / is a Planck 
function. 

The fluctuations in the photon distribution function inherit an energy 
dependence edf/de from the source terms in the Boltzmann equation 
Separating out the background contribution to /, and its energy dependence, 
we can write 

f{r,, x, e, e) = /(e) [1 - 6(r?, x, e)d In f/d In e] , (11) 

so that the CMB spectrum is Planckian but with a direction-dependent tem- 
perature AT/T — 9. Using the Lorentz invariance of /, it is not difficult to 
show that the quadrupole and higher moments of are gauge-invariant. If we 
now substitute for / in (|5J|, we find the Boltzmann equation for O: 

d ( ® + ^ + e • V(9 + i>) = -an e a T e + —an c a T I 6(e')[(l + (e • e'f] de 1 
or, 1d7t J 

+ an e crT e ■ v b + <t> + tp ■ (12) 



Anisotropics in the Cosmic Microwave Background 7 



The formal solution of this equation is an integral along the line of sight 
n = — e, 

{Q(n) + ^} R = e- T {e(h) + ij)} E + [ e- T Sd77, (13) 

J E 

where R is the reception event, E is the emission event, and t = J an e aT drj is 
the optical depth back from R. The source term S is given by the right-hand 
side of (|12|l . but with replaced by —ip in the first term. 

We gain useful insight into the physics of anisotropy formation by approxi- 
mating the last scattering surface as sharp (which is harmless on large angular 
scales), and ignoring the quadrupole CMB anisotropy at last scattering. In this 
case lfH5|) reduces to 

f R ■ ■ 

[0(n) + i/j] r = Qq\ e + ip\ E -n- v h \ B + (tp + ^dij, (14) 

JE 

where ©o is the isotropic part of 0, and is proportional to the fluctuation in 
the photon energy density. The various terms in this equation have a simple 
physical interpretation. The temperature received along direction n is the 
isotropic temperature of the CMB at the last scattering event on the line 
of sight, 0o, corrected for the gravitational redshift due to the difference in 
potential between E and R, and the Doppler shift e • Vh\ E resulting from 
scattering off moving electrons. Finally, there is an additional gravitational 
redshift contribution arising from evolution of the gravitational potentials 3 . 



Machinery for an Accurate Calculation 

An accurate calculation of the CMB anisotropy on all scales where linear 
perturbation theory is valid requires a full numerical solution of the Boltz- 
mann equation. The starting point is to expand 0(0, x, e) in appropriate 
basis functions. For scalar perturbations, these are the contraction of the (ir- 
reducible) trace-free tensor products e^ 1 . . . e 1 ^ (the angle brackets denoting 
the trace-free part) with trace-free (spatial) tensors derived from derivatives 
of scalars |28M31ll33) . Fourier expanding the scalar functions, we end up form- 
ing contractions between e^ 1 . . . e 1 ^ and fc/jj . . . where k is the wavevector. 
These contractions reduce to Legendre polynomials of k-e, and so the normal- 
mode expansion of for scalar perturbations takes the form 

It is straightforward to show that the implied azimuthal symmetry about the 
wavevector is consistent with the Boltzmann equation i|12[l . Inserting the ex- 
pansion of into this equation gives the Boltzmann hierarchy for the moments 
; : 
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((5/o - 1)6/ - Snv h + ^9 2 



+ 5i^ + 8nk^, (16) 

where v\, — J ikv\ } (k)c lk x d 3 fc/(27r) 3 / 2 , and cf> and ■0 are the Fourier trans- 
forms of the potentials. This system of ordinary differential equations can 
be integrated directly with the linearised Einstein equations for the metric 
perturbations, and the fluid equations governing perturbations in the other 
matter components, as in the publically-available COSMICS code [33] • Careful 
treatment of the truncation of the hierarchy is necessary to avoid unphysical 
reflection of power back down through the moments. 

A faster way to solve the Boltzmann equation numerically is to use the 
line-of-sight solution 113|) . as in the widely-used CMBFAST code j3J)] and 
its parallelised derivative CAMB [37]. Inserting the expansion (|15|l gives the 
integral solution to the hierarchy 

rvo r 

9/ 1„ = {21 + 1) I d V e- T (</> + (kA v ) - f(6 + i/>)ji (kA v ) 



,(17) 



where A-q = r) — 77, j\ is a spherical Bessel function, and primes denote 
derivatives with respect to the argument. Using the integral solution, it is 
only necessary to evolve the Boltzmann hierarchy to modest I to compute ac- 
curately the source terms that appear in the integrand. The integral approach 
is thus significantly faster than a direct solution of the hierarchy. 

The spherical multipoles a/ m of the temperature anisotropy can be ex- 
tracted from i|15|l as 

W / ^4We lfe - ■ (18) 

Statistical homogeneity and isotropy imply that the equal-time correlator 

(9/(77, k)Qt( v , fc')> = ^r '^ k ) s ( k ~ k ') . (1 9 ) 

so forming the correlation (aima*i m r) gives the power spectrum 

4tt 



Cr 



(21 +1) 



— Jef(k)dlnk. (20) 



If we consider (pure) perturbation modes characterised by a single indepen- 
dent stochastic amplitude per Fourier mode (such as the comoving curvature 
for the adiabatic mode; see Sect. 12 .4|) . the power 9^(fc) is proportional to 
the power spectrum of that amplitude. The spherical Bessel functions in ltT7|) 
peak sharply at kAr) = I for large I, so that multipoles I are mainly probing 
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spatial structure with wavenumber k ~ l/Arj at last scattering. The oscilla- 
tory tails of the Bessel functions mean that some power from a given k does 
also enter larger scale anisotropies. Physically, this arises from Fourier modes 
that are not aligned with their wavevector perpendicular to the line of sight. 
As we discuss in the next section, the tightly-coupled system of photons and 
baryons undergoes acoustic oscillations prior to recombination on scales in- 
side the sound horizon. For the pure perturbation modes, all modes with a 
given wavenumber reach the maxima or minima of their oscillation at the 
same time, irrespective of the direction of fc, and so we expect modulation in 
the C;s on sub-degree scales. The first three of these acoustic peaks have now 
been measured definitively; see Fig. ^ 

2.4 Photon Baryon Dynamics 

Prior to recombination, the mean free path of CMB photons is <~ 4.9 x 
10 4 (ilb/i 2 ) -1 (1 + z)~ 2 Mpc. On comoving scales below this length the pho- 
tons and baryons behave as a tightly-coupled fluid, with the CMB almost 
isotropic in the baryon frame. In this limit, only the I = and I — 1 moments 
of the distribution function are significant. 

The stress-energy tensor of the photons is given in terms of the distribution 
function by 

T»» = a' 2 J f( V , x, e, e)p»p v e dede , (21) 

so that the Fourier modes of the fractional over-density of the photons are 
6ry = 49o and the photon (bulk) velocity w 7 = — Oi. The anisotropic stress is 
proportional to 02. In terms of these variables, the first two moment equations 
of the Boltzmann hierarchy become 

, (22) 

t(v 7 - «b) . (23) 

Here, the derivative of the optical depth f = — an e CTx (and so is negative). 
The momentum exchange between the photons and baryons due to the drag 
term in ill'MI) gives rise to a similar term in the Euler equation for the baryons: 

Vb+Hv b + ki() = R~ 1 T(v b -v 7 ) , (24) 

where we have ignored baryon pressure. The ratio of the baryon energy density 
to the photon enthalpy is R = 3pb/4/3 7 and is proportional to the scale factor 
a, and TL = a/ a is the conformal Hubble parameter. 

In the tightly-coupled limit | -r 1 1 <C and Hr 1 . In this limit, we can 
treat the ratios of the mean-free path to the wavelength and the Hubble time 
as small perturbative parameters. Equations l|23|l and (|24|) then imply that 



40 



1 2 

— kSy — —k<d2 + kip 
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v ~t = to first order in the small quantities k/\f\ and H/\f\. Comparing the 
continuity equation for the baryons, 

5 h - kv h - 30 = , (25) 

with that for the photons, we see that 5 7 = 45b/3, so the evolution of the 
photon-baryon fluid is adiabatic, preserving the local ratio of the number 
densities of photons to baryons. Combining ill' Ml) and (|24|l to eliminate the 
scattering terms, and then using v 7 = v\>, we find the evolution of the photon 
velocity to leading order in tight coupling: 

7)7 + iTR 71 " 1 + 4(1 + R) kSl + k ^ = °- ( 26 ) 

The I > 1 moments of the photon distribution function arise from the bal- 
ance between isotropisation by scattering and their generation by photons free 
streaming over a mean free path; these moments are suppressed by factors 
(fc/|f|) i_1 . In particular, during tight coupling 9 2 « (20/27)fcf _1 -y 7 ignoring 
polarization. (The factor 20/27 rises to 8/9 if we correct for polarization |38|.) 

Combining Ij26(l with the photon continuity equation 1221) shows that the 
tightly-coupled dynamics of <5 7 is that of a damped, simple-harmonic oscillator 
driven by gravity |39| : 



S nR i 1 ,2r A 1 ^HR ; 4, 2 



The damping term arises from the redshifting of the baryon momentum in an 
expanding universe, while photon pressure provides the restoring force which 
is weakly suppressed by the additional inertia of the baryons. The WKB 
solutions to the homogeneous equation are 

S y = (1 + i?)~ 1/4 cosfcr s , and <5 7 = (1 + R)- 1/4 cos kr s , (28) 



where the sound horizon r s = JJ' Arf / yi{T+R). Note also that for static 
potentials, and ignoring the variation of R with time, the mid-point of the 
oscillation of <5 7 is shifted to —4(1 + R)ij). The dependence of this shift on the 
baryon density produces a baryon-dependent modulation of the height of the 
acoustic peak in the temperature anisotropy power spectrum; see Section [3] 

The driving term in l|27|) depends on the evolution of the gravitational 
potentials. If we ignore anisotropic stress, <f) and ip are equal, and their Fourier 
modes evolve as 



3H(1+tU- 



2H + I 1 + 3^ I H 2 



'-k' 
P 

^(Sp-Vjp) (L>!)) 



in a flat universe, which follows from the perturbed Einstein field equations. 
Here, p and p are the total density and pressure in the background model, 
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Fig. 2. Evolution of the potential <f> in adiabatic and CDM-isocurvature models 
for wavenumbers k — 0.001, 0.01 and 0.1 Mpc -1 (top to bottom respectively in 
matter domination). The conformal time at matter-radiation equality rj cq and last 
scattering jy* are marked by arrows. 



bp and bp are the Fourier modes of their perturbations, and n = 8nG. The 
source term is gauge- invariant; it vanishes for mixtures of barotropic fluids 
[pi = Pi(pi)] with 5pi/{pi + pi) the same for all components. For adiabatic 
perturbations, this latter condition holds initially and is preserved on super- 
Hubble scales. It is also preserved in the tightly-coupled photon-baryon fluid 
as we saw above. For adiabatic perturbations, the potential is constant on 
scales larger than the sound horizon when p/p is constant, but decays during 
transitions in the equation of state, such as from matter to radiation dom- 
ination. Above the sound horizon in flat models, it can be shown that the 
quantity 

H=-(j>- 2—2-; (30) 

na z (p + pj 

is conserved even through such transitions. The perturbation to the intrinsic 
curvature of comoving hypersurfaces (i.e. those perpendicular to the the four- 
velocity of observers who see no momentum density) is given in terms of 1Z as 
4(fc 2 / 'a 2 )TZ. Using the constancy of 1Z on large scales, the potential falls by a 
factor of 9/10 during the transition from radiation to matter domination. The 
evolution of the potential is illustrated in Fig. [21 in a flat ACDM model with 
parameters given in Sect. ^ The potential oscillates inside the sound hori- 
zon during radiation domination since the photons, which are the dominant 
component at that time, undergo acoustic oscillations on such scales. 

The behaviour of the potentials for isocurvature perturbations is quite dif- 
ferent on large scales during radiation domination , since the source term 
in 129|) is then significant. In isocurvature fluctuations, the initial perturba- 
tions in the energy densities of the various components compensate each other 
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Fig. 3. Evolution of the combination <5 7 /4 + ip (top left) and the photon velocity 
v 1 (bottom left) which determine the temperature anisotropies produced at last 
scattering (denoted by the arrow at 77*). Three modes are shown with wavenumbers 
k — 0.001, 0.1 and 0.2 Mpc -1 , and the initial conditions are adiabatic. The fluctu- 
ations at the time of last scattering are shown as a function of linear scale in the 
right-hand plot. 



in such a way that the comoving curvature 1Z — 0. Figure [21 shows the evolu- 
tion of CDM-isocurvature modes, in which there is initially a large fractional 
perturbation in the dark matter density, with a small compensating fractional 
perturbation in the radiation. (The full set of possibilities for regular isocur- 
vature modes are discussed in 03-) On large scales in radiation domination 
the potential grows as a, the scale factor. 

Adiabatic Fluctuations 

For adiabatic fluctuations, the photons are initially perturbed by <5 7 (0) = 
— 2 , 0(O) = 47\L(0)/3, i.e. they are over-dense in potential wells, and their ve- 
locity vanishes v 7 (0) = 0. If we consider super-Hubble scales at last scatter- 
ing, there has been insufficient time for v 7 to grow by gravitational infall and 
the action of pressure gradients and it remains small. The photon continuity 
equation lll'L'l) then implies that 5 1 — 4(f) remains constant, and the decay of 
4> through the matter-radiation transition leaves (<5 7 /4 + ip)^*) ~ ( / ) ( 7 7*)/3 = 
— 37?.(0)/5 on large scales (k < 3 x 10 _3 Mpc _1 ) at last scattering. The com- 
bination <5 7 /4 + ip = Oo + ip is the dominant contribution to the large-scale 
temperature anisotropies produces at last scattering; see (|14H . The evolution 
of the photon density and velocity perturbations for adiabatic initial condi- 
tions are show in Fig. |3 along with the scale dependence of the fluctuations 
at last scattering. The plateau in (<5 7 /4 + ip)(r]*) on large scales ensures that 
a scale-invariant spectrum of curvature perturbations translates into a scale- 
invariant spectrum of temperature anisotropies, 1(1 + l)Ci = constant, for 
small I. 

On scales below the sound horizon at last scattering, the photon-baryon 
fluid has had time to undergo acoustic oscillation. The form of the photon ini- 
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Fig. 4. Contribution of the various terms in 1141 to the temperature-anisotropy 
power spectrum from adiabatic initial conditions. At high I, the contributions are 
(from top to bottom): total power; S-,/4 + ip (denoted SW for Sachs-Wolfe [3]); 
Doppler effect from Vb\ and the integrated Sachs- Wolfe effect (1SW) coming from 
evolution of the potential along the line of sight. 

tial condition, and the observation that the driving term in (|27|) mimics the 
cosine WKB solution of the homogeneous equation (see Fig. |2J, set the oscil- 
lation mostly in the cos kr s mode. The midpoint of the oscillation is roughly 
at <5 7 /4 = —(1 + R)ijj. This behaviour is illustrated in Fig. Modes with 
kr s (j]*) = 7T have undergone half an oscillation at last scattering, and are max- 
imally compressed. The large value of Oo + tp at this particular scale gives rise 
to the first acoustic peak in Fig.^ now measured to be at I — 220.1 ±0.8 |4*2| . 
The subsequent extrema of the acoustic oscillation at fcr s (?7*) = nix give rise 
to the further acoustic peaks. The angular spacing of the peaks is almost 
constant and is set by the sound horizon at last scattering and the angular 
diameter distance to last scattering. The acoustic part of the anisotropy spec- 
trum thus encodes a wealth of information on the cosmological parameters; 
see Sect. |3 The photon velocity v 1 oscillates as sinfcr s , so the Doppler term 
in (|14fl tends to fill in power between the acoustic peaks. The relative phase 
of the oscillation of the photon velocity has important implications for the 
polarization properties of the CMB as discussed in Sect.@J The contributions 
of the various terms in (|14|) to the temperature-anisotropy power spectrum 
are shown in Fig. 0] for adiabatic perturbations. 
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Isocurvature Fluctuations 

For the CDM-isocurvature mode 1 the photons are initially unperturbed, as 
is the geometry: <$ 7 (0) = = 0(0) and v 1 = 0. On large scales <5 7 /4 = <f) is 
preserved, so the growth in <f> during radiation domination is matched by a 
growth in 8 1 and the photons are under-dense in potential wells. It follows that 
at last scattering (<5 7 /4 + ip)(r}*) ~ 20(?7*) for k < 3 x 10~ 3 Mpc -1 . Note that 
the redshift climbing out of a potential well enhances the intrinsic temperature 
fluctuation due to the photon under-density there. The evolution of the photon 
fluctuations for isocurvature initial conditions are shown in Fig. 

The evolution of the potential for isocurvature modes makes the driving 
term in (|27|l mimic the sine solution of the homogeneous equation, and so 5 7 
follows suit oscillating as sin ~ kr s about the equilibrium point —4(1 + R)ip. 
The acoustic peaks are at fcr s (?7*) ~ rt7r/2, and the photons are under-dense 
in the potential wells for the odd-n peaks, while over-dense in the even n. 
The various contributions to the temperature-anisotropy power spectrum for 
isocurvature initial conditions are shown in Fig. El The different peak po- 
sitions for isocurvature initial conditions allow the CMB to constrain their 
relative contribution to the total fluctuations. Current constraints are rather 
dependent on whether one allows for correlations between the adiabatic and 
isocurvature modes (as are generic in the multi-field inflation models that 
might have generated the initial conditions), and the extent to which addi- 
tional cosmological constraints are employed; see 03] for a recent analysis 
allowing for the most general correlations but a single power-law spectrum. 

1 It is also possible to have the dominant fractional fluctuation in the baryon density 
rather than the cold dark matter. However, this mode is nearly indistinguishable 
from the CDM mode since, in the absence of baryon pressure, they differ only by 
a constant mode in which the radiation and the geometry remain unperturbed, 
but the CDM and baryon densities have compensating density fluctuations |43|. 
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Fig. 6. As Fig.^Jbut for CDM-isocurvature initial conditions. The initial spectrum 
of entropy perturbations is scale-invariant. 

Beyond Tight-Coupling 

On small scales it is necessary to go beyond tight-coupling of the photon- 
baryon system since the photon diffusion length can become comparable to the 
wavelength of the fluctuations. Photons that have had sufficient time to diffuse 
of the order of a wavelength can leak out of over-densities, thus damping the 
acoustic oscillations and generating anisotropy 0^1 • A rough estimate of the 
comoving scale below which diffusion is important is the square root of the 
geometric mean of the particle horizon (or conformal age) and the mean-free 
path of the photons, i.e. \J rj/\f\. Converting this to a comoving wavenumber 
defines the damping scale 

fc D 2 ~ 0.3(an^ 2 ) _1/2 (^b/i 2 ) _1 (a/a») 5/2 Mpc 2 (31) 

when the scale factor is a. Here, a» is the scale factor at last scattering, and 
the expression is valid well after matter-radiation equality but well before 
recombination. The effect of diffusion is to damp the photon (and baryon) 
oscillations exponentially by the time of last scattering on comoving scales 
smaller than ~ 3 Mpc. The resulting damping effect on the temperature power 
spectrum has now been measured by several experiments (201 HH 1231 • 

To describe diffusion damping more quantitatively, we consider scales that 
were already sub-Hubble during radiation domination. The gravitational po- 
tentials will then have been suppressed during their oscillatory phase when the 
photons (which are undergoing acoustic oscillations themselves) dominated 
the energy density, and so we can ignore gravitational effects. Furthermore, 
the dynamical timescale of the acoustic oscillations is then short compared to 
the expansion time and we can ignore the effects of expansion. In this limit, 
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the Euler equations for the photons and the baryons can be iterated to give 
the relative velocity between the photons and baryons to first order in k/\f\: 

(1 +iT 1 )0 7 - Vb) = Jfcf- 1 ^ . (32) 

Using momentum conservation for the total photon-baryon system gives 

1 2 

i) 7 + Rv h + -k5j - -k@ 2 = , (33) 

which can be combined with the derivative of (|32|l to give a new Euler equation 
for the photons correct to first order in tight coupling: 

(1 + R)v 7 « -hs, + -Jt—^kT-H, + ^k 2 r~\ . (34) 

Here, we have used 6 2 « 8fcf _1 Uy/9 which includes the correction due to po- 
larization. In the limit of perfect coupling, l|34l) reduces to l|2t)|) on small scales. 
The continuity equation for the photons, <L = 4fci> 7 /3 (+40), shows that the 
last two terms on the right of l|34l) are drag terms, and on differentiating gives 

kH- 1 {16 R 2 \ • fc 2 . n 

s ^WTr1{t 5 + TTr) 5 '' + WTr) s ^ - (35) 

The WKB solution is 

is the damping scale. 

The finite mean-free path of CMB photons around last scattering has an 
additional effect on the temperature anisotropics. The visibility function — re T 
has a finite width ~ 80 Mpc and so along a given line of sight photons will 
be last scattered over this interval. Averaging over scattering events will tend 
to wash out the anisotropy from wavelengths short compared to the width of 
the visibility function. This effect is described mathematically by integrating 
the oscillations in the spherical Bessel functions in (|17(l against the product 
of the visibility function and the (damped) perturbations. 

Boltzmann codes such as CMBFAST and CAMB [371 use the tight- 
coupling approximation at early times to avoid the numerical problems asso- 
ciated with integrating the stiff Euler equations in their original forms 1231) 
and 

2.5 Other Features of the Temperature- Anisotropy Power 
Spectrum 

We end this section on the fundamentals of the physics of CMB temperature 
anisotropics by reviewing three additional effects that contribute to the linear 
anisotropics. 
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Integrated Sachs Wolfe Effect 

The integrated Sachs- Wolfe (ISW) effect is described by the last term on the 
right of 114(1. It is an additional source of anisotropy due to the temporal 
variation of the gravitational potentials along the line of sight: if a potential 
well deepens as a CMB photon crosses it then the blueshift due to infall will 
be smaller than redshift from climbing out of the (now deeper) well. (The 
combination <f> + ip has a direct geometric interpretation as the potential for 
the electric part of the Weyl tensor The ISW receives contributions 

from late times as the potentials decay during dark-energy domination, and 
at early times around last scattering due to the finite time since matter- 
radiation equality. 

The late-time effect contributes mainly on large angular scales since there 
is little power in the potentials at late times on scales that entered the Hubble 
radius during radiation domination. The late ISW effect is the only way to 
probe late-time structure growth (and hence e.g. distinguish between different 
dark-energy models) with linear CMB anisotropies, but this is hampered by 
cosmic variance on large angular scales. The late ISW effect produces corre- 
lations between the large-scale temperature fluctuations and other tracers of 
the potential in the local universe, and with the advent of the WMAP data 
these have now been tentatively detected 03011021; see also Sect. El 

In adiabatic models the early-time ISW effect adds coherently with the 
contribution <5 7 /4 + ip to the anisotropies near the first peak, boosting this 
peak significantly see Fig. 01 The reason is that the linear scales that 
contribute here are maximally compressed with <5 7 /4 + ip ~ —ip/2 which has 
the same sign as <f> for decaying (f>. 

Reionization 

Once structure formation had proceeded to produce the first sources of ultra- 
violet photons, the universe began to reionize. The resulting free electron 
density could then re-scatter CMB photons, and this tended to isotropise the 
CMB by averaging the anisotropies from many lines of sight at the scattering 
event. Approximating the bi-modal visibility function as two delta functions, 
one at last scattering 2 77* and one at reionization r) Te , if the optical depth 
through reionization is r ro , the temperature fluctuation at x = at rj is 

[9(n) + V]r, » (1 - e~ rrc )(9 + ip - n ■ v h )[-h(r) - rj IC ),Vrc) 

+ c~ Trc (6 + tp-h- v h )[-n(rj - V*),V*} ■ (37) 

Here, we have used (fHfl) . neglected the ISW effect, and approximated the scat- 
tering as isotropic. The first term on the right describes the effect of blending 

2 We continue to refer to the last scattering event around recombination as last 
scattering, even in the presence of re-scattering at reionization. 



18 Anthony Challinor 



the anisotropies from different lines of sight (to give 80) and the generation of 
new anisotropies by re-scattering off moving electrons at reionization; the sec- 
ond term is simply the temperature anisotropy that would be observed with 
no reionization, weighted by the fraction of photons that do not re-scatter. 
Since 80 + ip at the re-scattering event is the average of 8o + 1/> — n -v^ on the 
electron's last scattering surface, on large scales k(rj le — 77*) <C 1 it reduces to 
80 + ip at [— n(rja — ?7* ),?/*], while on small scales it vanishes. It follows that 
for scales that are super-horizon at reionization, the observed temperature 
anisotropy becomes 

8(n) -> 0(n) - (1 -e- r ")n • Av h , (38) 

where Av^ is the difference between the electron velocity at the reionization 
event and the preceding last scattering event on the line of sight. On such 
scales the Doppler terms do not contribute significantly and the temperature 
anisotropy is unchanged. For scales that are sub-horizon at reionization, 

0(n) -» c- r "-'8(n) - (1 - c- r ")n • v h , (39) 

where the Doppler term is evaluated at reionization. In practice, the visibility 
function is not perfectly sharp at reionization and the integral through the fi- 
nite re-scattering distance tends to wash out the Doppler term since only plane 
waves with their wavevectors near the line of sight contribute significantly to 
n ■ v\,. Figure shows the resulting effect C; — > e _2rro C; on the anisotropy 
power spectrum on small scales. Recent results from WMAP \Ftj\ suggest an 
optical depth through reionization r re ~ 0.17. Such early reionization cannot 
have been an abrupt process since the implied redshift z rc ~ 15 is at odds with 
the detection of traces of smoothly-distributed neutral hydrogen at z ~ 6 via 
Gunn-Peterson troughs in the spectra of high-redshift quasars |5(JI 15 1| . 



Tensor Modes 



Tensor modes, describing gravitational waves, represent the transverse trace- 
free perturbations to the spatial metric: 



As 2 = a 2 (ri)[dr] 2 - (<% + Ti^dscMx 3 '] , 



(40) 



with h\ = and dM = 0. A convenient parameterisation of the photon 
four-momentum in this case is 



he* 



-h i e j 
2 3 



(41) 



where e 2 = 1 and e is a times the energy of the photon as seen by an observer 
at constant x. The components of e are the projections of the photon direction 
for this observer on an orthonormal spatial triad of vectors a~ 1 (9 i — hjdj/2). 
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Fig. 7. Effect of reionization on the temperature-anisotropy power spectrum. The 
spectra are (from top to bottom) for no reionization, r rc = 0.1 and 0.2. 



In the background e = da;/d?7 and is constant. The evolution of the comoving 
energy e in the perturbed universe is 



1 de 1 • 

e cby + 2 



+ -AyeV = , (42) 



and so the Boltzmann equation for 0(??, x, e) is 

^ + e • V6 = -an c <j T e + —an c a T [ 6(e')[(l + (e • e') 2 ] de' 
or] 16-k J 

- . (43) 

Neglecting the anisotropic nature of Thomson scattering, the solution of this 
equation is an integral along the unperturbed line of sight: 



1 



6(n) = / e- T h l3 n l n J drj . (44) 

The time derivative hij is the shear induced by the gravitational waves. This 
quadrupole perturbation to the expansion produces an anisotropic redshifting 
of the CMB photons and an associated temperature anisotropy. 

Figure [S] compares the power spectrum due to gravitational waves with 
that from scalar perturbations for a tensor-to-scalar ratio r = 1 corresponding 
to an energy scale of inflation 3.3 x 10 16 GeV. The constraints on gravitational 
waves from temperature anisotropies are not very constraining since their ef- 
fect is limited to large angular scales where cosmic variance from the dominant 
scalar perturbations is large. Gravitational waves damp as they oscillate inside 
the horizon, so the only significant anisotropies are from wavelengths that are 
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Fig. 8. The temperature-anisotropy power spectrum from scalar perturbations (den- 
sity perturbations; top) and tensor perturbations (gravity waves; bottom) for a 
tensor-to-scalar ratio r = 1. 

super-horizon at last scattering, corresponding to I ~ 60. The current 95-per 
cent upper limit on the tensor-to-scalar ratio is 0.68 20.. Fortunately, CMB 
polarization provides an alternative route to detecting the effect of gravita- 
tional waves on the CMB which is not limited by cosmic variance [5211^3; see 
also Sect. |U 

3 Cosmological Parameters and the CMB 

The simple, linear physics of CMB temperature anisotropics, reviewed in the 
previous section, means that the CMB depends sensitively on many of the 
key cosmological parameters. For this reason, CMB observations over the 
past decade have been a significant driving force in the quest for precision 
determinations of the cosmological parameters. It is not our intention here 
to give a detailed description of the constraints that have emerged from such 
analyses, e.g. JU], but rather to provide a brief description of how the key 
parameters affect the temperature-anisotropy power spectrum. More details 
can be found in the seminal papers on this subject, e.g. EHI 02] an d 
references therein. 

3.1 Matter and Baryons 

The curvature of the universe and the properties of the dark energy are largely 
irrelevant for the pre-recombination physics of the acoustic oscillations. Their 
main contribution is felt geometrically through the angular diameter distance 
to last scattering, Da, which controls the projection of linear scales there to 
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Fig. 9. Dependence of the temperature-anisotropy power spectrum on the physical 
density in baryons (left) and all non-relativistic matter (right). From top to bottom 
at the first peak, the baryon densities vary linearly in the range Q. h h 2 = 0.06-0.005 
(left) and the matter densities in f2 m /i 2 = 0.05-0.5 (right). The initial conditions 
are adiabatic. 

angular scales on the sky today. In contrast, those parameters that determine 
the energy content of the universe before recombination, such as the physical 
densities in (non-relativistic) matter 17 m /i 2 , and radiation fl T h 2 (determined by 
the CMB temperature and the physics of neutrinos), play an important role in 
acoustic physics by determining the expansion rate and hence the behaviour of 
the perturbations. In addition, the physical density in baryons, f2b/i 2 , affects 
the acoustic oscillations through baryon inertia and the dependence of the 
photon mean-free path on the electron density. The effect of variations in the 
physical densities of the matter and baryon densities on the anisotropy power 
spectrum is illustrated in Fig. for adiabatic initial conditions. 

The linear scales at last scattering that have reached extrema of their 
oscillation are determined by the initial conditions (i.e. adiabatic or isocurva- 
ture) and the sound horizon r s (r)*). Increasing the baryon density holding the 
total matter density fixed reduces the sound speed while preserving the ex- 
pansion rate (and moves last scattering to slightly earlier times) . The effect is 
to reduce the sound horizon at last scattering and so the wavelength of those 
modes that are at extrema of their oscillation, and hence push the acoustic 
peaks to smaller scales. This effect could be confused with a change in the 
angular diameter distance Da, but fortunately baryons have another distin- 
guishing effect. Their inertia shifts the zero point of the acoustic oscillations 
to ~ —(1 + R)ip, and enhances the amplitude of the oscillations. In adiabatic 
models for modes that enter the sound horizon in matter domination, <5 7 /4 
starts out at —2ip/3, and so the amplitude of the oscillation is — + 3i?)/3. 
The combination of these two effects is to enhance the amplitude of 6o + t/> 
at maximal compression by a factor of 1 + 6R over that at minimal compres- 
sion. The effect on the power spectrum is to enhance the amplitude of the 
1st, 3rd etc. peaks for adiabatic initial conditions, and the 2nd, 4th etc. for 
isocurvature. Current CMB data gives fl^h 2 = 0.023 ± 0.001 for power-law 
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Fig. 10. Dependence of the temperature- anisotropy power spectrum on the curva- 
ture Q.k (left) and cosmological constant £7a (right) in adiabatic models. In both 
cases, the physical densities in baryons and matter were held constant, thus preserv- 
ing the conditions on the last scattering surface. The curvature varies (left to right) 
in the range -0.15-0.15 and the cosmological constant in the range 0.9-0.0. 



ACDM models , beautifully consistent with determinations from big bang 
nucleosynthesis. Other effects of baryons are felt in the damping tail of the 
power spectrum since increasing the baryon density tends to inhibit diffusion 
giving less damping at a given scale. 

The effect of increasing the physical matter density £l m h 2 at fixed Q.\ ) h 2 is 
also two-fold (see Fig.|5J): (i) a shift of the peak positions to larger scales due 
to the increase in Z?a! & n d (ii) a scale-dependent reduction in peak height in 
adiabatic models. Adiabatic modes that enter the sound horizon during radi- 
ation domination see the potentials decay as the photon density rises to reach 
maximal compression. This decay tends to drive the oscillation, increasing the 
oscillation amplitude. Raising fl m h 2 brings matter-radiation equality to ear- 
lier times, and reduces the efficiency of the gravitational driving effect for the 
low-order peaks. Current CMB data gives fl m h 2 = 0.13 ± 0.01 for adiabatic, 
power-law ACDM models ■ 

3.2 Curvature, Dark Energy and Degeneracies 

The main effect of curvature and dark energy on the linear CMB anisotropies 
is through the angular diameter distance and the late-time integrated Sachs- 
Wolfe effect; see Fig. ^1 for the case of adiabatic fluctuations in cosmological- 
constant models. The ISW contribution is limited to large scales where cosmic 
variance severely limits the precision of power spectrum estimates. There is 
an additional small effect due to quantisation of the allowed spatial modes 
in closed models (e.g. but this is also confined to large scales (i.e. near 

the angular projection of the curvature scale). Most of the information that 
the CMB encodes on curvature and dark energy is thus locked in the angular 
diameter distance to last scattering, Da- 
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Fig. 11. The geometric degeneracy. A scale-invariant adiabatic ACDM model with 
Q. h h 2 = 0.024, n m /i 2 = 0.14 and fi A = 0.73 and Q, K = (close to the WMAP 
best-fit values |12p produces an almost identical spectrum to a closed model Qk = 
—0.288 with vanishing cosmological constant. However, the Hubble constants are 
very different - h = 0.72 in the flat model and 0.33 in the closed model - and so 
the latter is easily ruled out by extern al constrain ts. The shaded region shows the 
lcr cosmic variance errors AC; jCi = y/2/(2l + 1) on the power spectrum. 



With the physical densities Q^ti 1 and £l m h 2 fixed by the acoustic part of 
the anisotropy spectrum, Da can be considered a function of VLk and the 
history of the energy density of the dark energy (often modelled through 
its current density and a constant equation of state). In cosmological con- 
stant models Da is particularly sensitive to the curvature: the 95-per cent 
interval from WMAP alone (with the weak prior Hq > 50 kms -1 Mpc -1 ) is 
—0.08 < £Ik < 0.02, so the universe is close to being spatially flat. The fact 
that the impact of curvature and the properties of the dark energy on the 
CMB is mainly through a single number Da leads to a geometrical degener- 
acy in parameter estimation |56) . as illustrated in Fig. 1111 Fortunately, this is 
easily broken by including other, complementary cosmological datasets. The 
constraint on curvature from WMAP improves considerably when supernovae 
measurements [37| or the measurement of Hq from the Hubble Space 
Telescope Key Project are included. Other examples of near-perfect de- 
generacies for the temperature anisotropies include the addition of gravity 
waves and a reduction in the amplitude of the initial fluctuations mimicing 
the effect of reionization. This degeneracy is broken very effectively by the 
polarization of the CMB. 
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4 CMB Polarization 

The growth in the mean-free path of the CMB photons during recombination 
allowed anisotropies to start to develop. Subsequent scattering of the radia- 
tion generated (partial) linear polarization from the quadrupole anisotropy. 
This linear polarization signal is expected to have an r.m.s. ~ 5 /iK, and, for 
scalar perturbations, to peak around multipoles I ~ 1000 corresponding to the 
angle subtended by the mean-free path around last scattering. The detection 
of CMB polarization was first announced in 2002 by the Degree Angular Scale 
Interferometer (DASI) team WMAP has also detected the polarization 
indirectly through its correlation with the temperature anisotropies ^Hj- A 
direct measurement of the polarization power from two-years of WMAP data 
is expected shortly. Polarization is only generated by scattering, and so is a 
sensitive probe of conditions at recombination. In addition, large-angle polar- 
ization was generated by subsequent re-scattering as the universe reionized, 
providing a unique probe of the ionization history at high redshift. 

4.1 Polarization Observables 

Polarization is conveniently described in terms of Stokes parameters I, Q, 
U and V, where / is the total intensity discussed at length in the previous 
section. The parameter V describes circular polarization and is expected to 
be zero for the CMB since it is not generated by Thomson scattering. The 
remaining parameters Q and U describe linear polarization. They are the 
components of the trace-free, (zero-lag) correlation tensor of the electric field 
in the radiation, so that for a quasi-monochromatic plane wave propagating 
along the z direction 




where the angle brackets represent an average on timescales long compared 
to the period of the wave. For diffuse radiation we define the polarization 
brightness tensor V a b(n) to have components given by ijlSl for plane waves 
within a bundle around the line of sight n and around the specified frequency. 
The polarization tensor is transverse to the line of sight, and, since it inherits 
its frequency dependence from the the quadrupole of the total intensity, has 
a spectrum given by the derivative of the Planck function (see equation II ljl . 

The polarization tensor can be decomposed uniquely on the sphere into 
an electric (or gradient) part and a magnetic (or curl) part [521 153) : 

Vab = V (a V fc) P B - e c {a V b) V c P B , (46) 

where angle brackets denote the symmetric, trace-free part, V a is the covariant 
derivative on the sphere, and e & is the alternating tensor. The divergence 
V a 7 , a 6 is a pure gradient if the magnetic part Pb = 0, and a curl if the electric 
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Fig. 12. Polarization patterns for a pure-electric mode (left) and pure-magnetic 
mode (right) on a small patch of the sky for potentials that are locally Fourier 
modes. The shading denotes the amplitude of the potential. For the electric pattern 
the polarization is aligned with or perpendicular to the Fourier wavevector depending 
on the sign of the potential; for the magnetic pattern the polarization is at 45 degrees. 



part Pe = 0. The potential Pe is a scalar under parity, but Pb is a pseudo- 
scalar. For a given potential P, the electric and magnetic patterns it generates 
(i.e. with P E = P and P B = P respectively) are related by locally rotating 
the polarization directions by 45 degrees. The polarization orientations on a 
small patch of the sky for potentials that are locally Fourier modes are shown 
in Fig. 1121 The potentials can be expanded in spherical harmonics (only the 
/ > 2 multipoles contribute to V a b) as 



P B (n) 



(The normalisation is conventional.) Under parity Ei m — > (— 1) Ei m but 
Bim — * — (— 1) Bi m . Assuming rotational and parity invariance, B is not 
correlated with E or the temperature anisotropics T, leaving four non- 
vanishing power spectra: Cj ', Cf , Cf and the cross-correlation Cf E , where 
e.g. (E lm T*J = Cj E . 





4.2 Physics of CMB Polarization 



For scalar perturbations, the quadrupole of the temperature anisotropies at 
leading order in tight coupling is 02 ~ kf~ 1 v 7 . Scattering of this quadrupole 
into the direction — n generates linear polarization parallel or perpendicular to 
the projection of the wavevector k onto the sky, i.e. V L j ~ 02[fc(ifcj\] TT , where 
TT denotes the transverse (to n), trace- free part. In a flat universe the po- 
larization tensor is conserved in the absence of scattering; for non-flat models 
this is still true if the components are defined on an appropriately-propagated 
basis (e.g. j2E])- For a single plane wave perturbation, the polarization on the 
sky is thus purely electric (see Fig. H2jl . For tensor perturbations, the polar- 
ization Vij ~ f _1 [/iy] TT since the tightly-coupled quadrupole is proportional 
to the shear hij. The gravitational wave defines additional directions on the 
sky when its shear is projected, and the polarization pattern is not purely 
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Fig. 13. Power spectra produced by adiabatic scalar perturbations (left) and tensor 
perturbations (right) for a tensor-to-scalar ratio r = 1. On large scales the spectra 
from scalar perturbations are (from top to bottom) Cf , Cj B and Of. For tensor 
perturbations, they are Cf , Cf E , Cf and Cf . 

electric. Thus density perturbations do not produce magnetic polarization in 
linear perturbation theory, while gravitational waves produce both electric 
and magnetic 

The polarization power spectra produced by scalar and tensor perturba- 
tions are compared in Fig.^3 The scalar Cf spectrum peaks around / ~ 1000 
since this corresponds to the projection of linear scales at last scattering for 
which diffusion generates a radiation quadrupole most efficiently. The polar- 
ization probes the photon bulk velocity at last scattering, and so Cf peaks 
at the troughs of Cj , while Cj E is zero at the peaks and troughs, and has 
its extrema in between. For adiabatic perturbations, the large-scale cross- 
correlation changes sign at I ~ 50, and, with the conventions adopted here 3 is 
positive between I = 50 and the first acoustic peak in Cj . Isocurvature modes 
produce a negative correlation from I = 2 to the first acoustic trough. 

Tensor modes produce similar power in electric and magnetic polarization. 
As gravitational waves damp inside the horizon, the polarization peaks just 
shortward of the horizon size at last scattering / ~ 100 despite these large 

3 The sign of Ei m for a given polarization field depends on the choice of conventions 
for the Stokes parameters and their decomposition into electric and magnetic 
multipoles. We follow |6U|. which produces the same sign of Cf E as H3, but note 
that the Boltzmann codes CMBFAST and CAMB |371 have the opposite sign. 
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scales being geometrically less efficient at transferring power to the quadrupole 
during a mean-free time than smaller scales. 

For both scalar and tensor perturbations, the polarization would be small 
on large scales were it not for reionization, since a significant quadrupole is 
only generated at last scattering when the mean-free path approaches the 
wavelength of the fluctuations. However, reionization does produce significant 
large-angle polarization (see Fig. ll3|) . The temperature quadrupole at last 
scattering peaks on linear scales with fc(?7 r e — r/») ~ 2, which then re-projects 
onto angular scales Z ^ 2(770 — rj r c)/(Vrc —f]*)- The position of the reionization 
feature is thus controlled by the epoch of reionization, and the height by the 
fraction of photons that scatter there i.e. r ro . The measurement of r re with 
large-angle polarization allows an accurate determination of the amplitude 
of scalar fluctuations from the temperature-anisotropy power spectrum. In 
addition, the fine details of the large-angle polarization power can in principle 
distinguish different ionization histories with the same optical depth, although 
this is hampered by the large cosmic variance at low I |62j . 

5 Highlights of Recent Results 

In this section we briefly review some of the highlights from recent observations 
of the CMB temperature and polarization anisotropies. Analysis of the former 
have entered a new phase with the release of the first year data from the 
WMAP satellite [H]; a further three years worth of data are expected from 
this mission. Detections of CMB polarization are still in their infancy, but 
here too we can expect significant progress from a number of experiments in 
the short term. 

5.1 Detection of CMB Polarization 

The first detection of polarization of the CMB was announced in September 
2002 X3] . The measurements were made with DASI, a compact interferomet- 
ric array operating at 30 GHz, deployed at the South Pole. The DASI team 
constrained the amplitude of the E and -B-mode spectra with assumed spec- 
tral shapes derived from a concordant ACDM model. They obtained a ~ 5-er 
detection of a non-zero amplitude for E with a central value perfectly consis- 
tent with that expected from the amplitude of the temperature anisotropies. 
DASI also detected the temperature-polarization cross-correlation at 95-per 
cent significance, but no evidence for B-mode polarization was found. The 
DASI results of a maximum-likelihood band-power estimation of the E and 
TE power spectra are given in Fig. El 

Measurements of Cj E were also provided in the first-year data release 
from WMAP, although polarization data itself was not released. These re- 
sults are also shown in Fig. [21 The existence of a cross-correlation between 
temperature and polarization on degree angular scales provides evidence for 
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Fig. 14. Current measurements (as of February 2004) of Cf E (top) and Cf (bot- 
tom). The points with 1-a errors are from the first one- year data release from 
WMAP |fv{) . The error boxes are the fiat band-power results from DASI centred 
on the maximum-likelihood band power and spanning the 68-per cent interval. The 
solid lines are the predicted power from the best-fit model to all the WMAP data. 

the existence of super-horizon fluctuations on the last scattering surface at 
recombination. This is more direct evidence for such fluctuations than from 
the large-scale temperature anisotropies alone, since the latter could have 
been generated gravitationally all along the line of sight. The sign of the 
cross-correlation and the phase of its acoustic peaks relative to those in the 
temperature-anisotropy spectrum is further strong evidence for adiabatic fluc- 
tuations. The one surprise in the WMAP measurement of Cf E is the be- 
haviour on large scales. A significant excess correlation over that expected if 
polarization were only generated at recombination is present on large scales 
(/ < 20). The implication is that reionization occurred early, 11 < z ve < 30, 
giving a significant optical depth for re-scattering: r re = 0.17 ±0.04 at 68-per 
cent confidence. As mentioned in Sect. 12.51 reionization at this epoch is ear- 
lier than that expected from observations of quasar absorption spectra and 
suggests a complex ionization history. 

5.2 Implications of Recent Results for Inflation 

The generic predictions from simple inflation models are that: (i) the uni- 
verse should be (very nearly) spatially flat; (ii) there should be a nearly 
scale-invariant spectrum of Gaussian, adiabatic density perturbations giving 
apparently-super-horizon fluctuations on the last scattering surface; and (iii) 
there should be a stochastic background of gravitational waves with a nearly 
scale-invariant (but necessarily not blue) spectrum. The amplitude of the lat- 



Anisotropies in the Cosmic Microwave Background 29 



ter is a direct measure of the Hubble rate during inflation, and hence, in 
slow-roll models, the energy scale of inflation. 

As discussed in Sect. 13.21 the measured positions of the acoustic peaks 
constrains the universe to be close to flat. The constraint improves further 
with the inclusion of other cosmological data. There is no evidence for isocur- 
vature modes in the CMB, although the current constraints are rather weak 
if general, correlated modes are allowed in the analysis |44| . Several of the 
cosmological parameters for the isocurvature models most favoured by CMB 
data are violently at odds with other probes, most notably the baryon den- 
sity which is pushed well above the value inferred from the abundances of the 
light-elements. There is also no evidence for primordial non-Gaussianity in 
the CMB (see e.g. |E2) 4 . 

Within flat ACDM models with a power-law spectrum of curvature fluc- 
tuations, the spectral index is constrained by the CMB to be close to scale 
invariant ^2] i although the inclusion of the latest data from small-scale exper- 
iments, such as CBI [70] and VSA |7T], tends to pull the best fit from WMAP 
towards redder power-law spectra: e.g. n s = 0.97l ; 3 at 68-per cent confi- 
dence combining WMAP and VSA [JTj. Slow-roll inflation predicts that the 
fluctuation spectrum should be close to a power law, with a run in the spectral 
index that is second order in slow roll: dn s /dlnfc ~ (n s — l) 2 . TheWMAP team 
reported weak evidence for a running spectral index by including small-scale 
data from galaxy redshift surveys and the Lyman-a forest, but modelling un- 
certainties in the latter have led many to question the reliability of this result 
(e.g. [Z2|)- New data from CBI and VSA now provide independent evidence for 
running in Hat ACDM models at the 2-cr level from the CMB alone. This re- 
flects the tension between the spectral index favoured by the low-/ CMB data 
(which is anomalously low for I < 10, favouring bluer spectra) and the high-/ 
data from the interferometers. The evidence for running is weakened consider- 
ably with the inclusion of external priors from large-scale structure data. The 
best-fit values for the run in n s obtained with the CMB alone are uncomfort- 
ably large for slow-roll inflation models, and give low power on small scales 
that is difficult to reconcile with the early reionization implied by the WMAP 
polarization data. However, a recent analysis |73| argues that the evidence for 
running depends crucially on the techniques employed to estimate the low-/ 
power from WMAP data, and that the running is strongly suppressed if exact 
likelihood techniques are adopted. A definitive answer on whether departures 
from power-law spectra are significant must probably await further data on 
both large and small scales. 

The final prediction of slow-roll inflation - the generation of nearly scale- 
invariant background of gravitational waves - is yet to be verified. The current 

4 The WMAP data does appear to harbour some statistically-significant depar- 
tures from rotational invariance 65 66 67 68 ST?). The origin of these effects, 
i.e. primordial or systematic due to instrument effects or imperfect foreground 
subtraction, is as yet unclear. 
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limits on the tensor-to-scalar ratio are only weak: [7T] report r < 0.68 at 95-per 
cent confidence from all CMB data in general, non-flat, adiabatic ACDM mod- 
els. Despite this, observations are beginning to place interesting constraints 
on specific models of inflation in the r-n s plane [^] [7S] . Already, large-field 
models with power-law potentials steeper than V oc 6 are ruled out due to 
their red scalar spectra and comparatively large tensor-to-scalar ratio. Future 
programmes targeting i?-mode polarization may ultimately be able to detect 
gravitational waves down to an inflationary energy scale of a few x 10 15 GeV. 
Such observations will sharpen constraints in the r-n s plane considerably, and 
should allow fine selection amongst the many proposed models of inflation. 

5.3 Detection of Late-Time Integrated Sachs- Wolfe Effect 

The late-time ISW effect arises from the decay of the gravitational potentials 
once the universe becomes dark-energy dominated, and so should produce 
large-angle (positive) correlations between the CMB temperature anisotropics 
and other tracers of the potential in the local universe. With the advent of 
the WMAP data, a number of groups have reported the detection of such a 
correlation. In WMAP data was cross-correlated with data on the hard 
X-ray background (which is dominated by emission from active galaxies) from 
the HEAO-I satellite, and the number density of radio sources from the NVSS 
catalogue. In each case a positive correlation was detected at significance 3cr 
and 2.5cr respectively. The correlation with NVSS has also been carried out 
independently by the WMAP team 0Hjj who also note that the observed 
positive correlation can be used to rule out the closed, A = model model 
that is a good fit to the CMB data in isolation (see Fig. ITT|) . Several groups 
have now also detected the cross-correlation on large scales between the CMB 
and optical galaxy surveys, e.g. 

6 Conclusion 

The linear anisotropics of the cosmic microwave background have been stud- 
ied theoretically for over three decades. The physics, which is now well under- 
stood, employs linearised radiative transfer, general relativity, and hydrody- 
namics to describe the propagation of CMB photons and the evolution of the 
fluid constituents in a perturbed Friedmann-Robertson- Walker universe. A 
number of bold predictions have emerged from this theoretical activity, most 
notably the existence of acoustic peaks in the anisotropy power spectrum due 
to oscillations in the photon-baryon plasma prior to recombination. Observers 
have risen to the challenge of verifying these predictions, and their detection 
is proceeding at a staggering rate. The large-scale Sachs- Wolfe effect, acoustic 
peak structure, damping tail, late-time integrated Sachs-Wolfe effect, polar- 
ization and reionization signature have all been detected, and the first three 
have been measured in considerable detail. Already, the size and scale of these 
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effects is allowing cosmological models to be constrained with unprecedented 
precision. The results are beautifully consistent with almost-scale-invariant 
adiabatic initial conditions evolving passively in a spatially flat, ACDM uni- 
verse. 

Much work still remains to be done to exploit fully the information con- 
tained in the CMB anisotropies. The Planck satellite, due for launch in 2007, 
should provide definitive mapping of the linear CMB anisotropies, and a 
cosmic- variance limited measurement of the power spectrum up to multipoles 
I ~ 2000. This dataset will be invaluable in assessing many of the issues hinted 
at in the first-year release of WMAP data, such as the apparent lack of power 
on large scales and possible violations of rotational (statistical) invariance. 
Prior to Planck, a number of ground-based programmes should shed further 
light on the issue of whether departures from a power- law primordial spectrum 
are required on cosmological scales, and the implications of this for slow-roll 
inflation. In addition, these small-scale observations will start to explore the 
rich science of secondary anisotropies, due to e.g. scattering in hot clusters |7H] 
or bulk flows modulated by variations in the electron density in the reionized 
universe [ZZIIZHjj and the weak lensing effect of large-scale structure |75] . 

Detections of CMB polarization are in their infancy, but we can expect 
rapid progress on this front too. Accurate measurements of the power spectra 
of -E-mode polarization, and its correlation with the temperature anisotropies, 
can be expected from a number of ground and balloon-borne experiments, as 
well as from Planck. The ultimate goal for CMB polarimetry is to detect the 
i3-mode signal predicted from gravitational waves. This would give a direct 
measure of the energy scale of inflation, and, when combined with measure- 
ments of the spectrum density perturbations, place tight constraints on the 
dynamics of inflation. Plans are already being made for a new generation of 
polarimeters with the large numbers of detectors and exquisite control of in- 
strument systematics needed to detect the gravity-wave signal if the energy 
scale of inflation is around 10 16 GeV. Ultimately, confusion due to imperfect 
subtraction of astrophysical foregrounds and the effects of weak lensing on the 
polarization limit will limit the energy scales that we can probe with CMB 
polarization; see (HOI and references therein. 
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